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Abstract 

We studied the dynamics of tlie proteome of influenza vims A/PR/8/34 (HINI) infected Madin-Darby canine kidney cells up 
to 12 hours post Infection by mass spectrometry based quantitative proteomics using the approach of stable isotope 
labeling by amino acids in cell culture (SILAC). We identified 1311 cell proteins and, apart from the proton channel M2, all 
major virus proteins. Based on their abundance two groups of virus proteins could be distinguished being In line with the 
function of the proteins In genesis and formation of new virions. Further, the data indicate a correlation between the 
amount of proteins synthesized and their previously determined copy number inside the viral particle. We employed 
bloinformatic approaches such as functional clustering, gene ontology, and pathway (KEGG) enrichment tests to uncover 
co-regulated cellular protein sets, assigned the Individual subsets to their biological function, and determined their 
Interrelation within the progression of viral Infection. For the first time we are able to describe dynamic changes of the 
cellular and, of note, the viral proteome in a time dependent manner simultaneously. Through cluster analysis, time 
dependent patterns of protein abundances revealed highly dynamic up- and/or down-regulation processes. Taken together 
our study provides strong evidence that virus infection has a major impact on the cell status at the protein level. 
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introduction 

The evolution of viruses is accompanied by an opposing 
evolution through constant interaction with their host. The high 
risk of infection by viruses that continually adapt strategies to 
overcome the cellular antiviral-defense is exemplified in the case of 
influenza A viruses by seasonal as well as pandemic outbreaks with 
serious consequences for the human population. For decades, 
enormous efforts are going on to understand the molecular details 
of viral replication itself, the complex interplay between viruses 
and hosts, and the corollary for the host cell to combat viral 
infection. Virus replication is strongly dependent on the biochem- 
ical, physiological, and physical status of the infected host cell. This 
is due to the involvement of various distinct cellular processes, 
several resources, and competing cellular requirements, such as 
biosynthesis and transport machinery during the infection process 
[!]■ 

The envelope of influenza A virus contains two major surface 
proteins, hemagglutinin (HA) and neuraminidase (NA), and - with 
a low abundance - the proton channel M2. The matrix protein 1 
(Ml) forms a layer beneath the viral membrane, enveloping eight 
different RNA segments. These segments are associated with the 



nucleoprotein (NP) as well as the three polymerase subunits (PA, 
PBl, PB2) forming viral ribonucleoprotein complexes (vRNPs). 
Upon binding to the host cell surface, influenza virus enters the 
cell via endocytic routes. After acidification of the endosome 
lumen a conformational change of HA triggers fusion of the 
envelope with the endosomal membrane releasing the segmented 
virus genome for transport into the nucleus. The genome is further 
encoding for two regulatory proteins: the nonstructural protein 1 
(NSl) which is expressed in the host cell but is not a component of 
the virion itself, and the nuclear export protein (NEP, synonymous 
NS2) which is represented in the virus particle in small quantities 
[2] . NS 1 suppresses transport of host mRNA from the nucleus, 
thus, favoring the export of viral mRNA, while NEP mediates the 
nuclear export of viral vRNPs [3-5] . 

The invasive hijacking of the host cell machinery by the virus is 
associated with directed influence on gene expression of host cell 
proteins [6] . Recently, the mimicry of cellular signal sequences by 
viral short linear motifs (SLiMs) was assigned to play a key role for 
viral hijacking of host transport and biosynthesis. So far, certain 
virus and host cell components interacting with each other have 
been identified through yeast-two-hybrid assays, computational 
approaches, and genome-wide RNA interference screens [7-12]. 
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These partners include RNA-binding proteins, transport proteins, 
transcription factors, and proteins of the intra-cellular signaling 
pathways (NFkB, apoptosis, MAPK, and WNT). However, a 
comprehensive view on the complex infection process and its 
consequences in particular for the proteome is still lacking. 

To unravel the molecular basis of virus-host-interaction on a 
systems level systematic approaches are required to conceive the 
whole replication process. Such deep investigation enables the 
identification of the viral Achilles heel as a potential target for 
effective antiviral strategies including respective drugs [13]. 
Former in vitro studies and in silico analyses addressing this issue 
found various host genes being essential for the influenza A virus 
infection cycle [8,14—16], However, it remained open how 
modification of gene expression translates into dynamics of the 
proteome and regulatory networks. Previously, quantitative 
proteomic analysis of influenza A virus (HlNl) infected A549 
cells [17-19] and primary human bronchial airway epithelial cells 
[19] have been performed by stable isotope labeling by amino 
acids in ceU culture (SILAC). Coombs et al. [17] found that about 
360 out of almost 4700 characterized cytosolic proteins were 
diflerentially regulated upon infection. Those are involved, for 
example, in cell immunity and antigen presentation, metabolism, 
signal transduction, and transcription pathways. This analysis was 
done 24 hrs post infection (p.i.), where significant titers of progeny 
viruses of HlNl infected A549 cells have been observed [20]. 
Hence, the study provides important information on the cell 
proteome in the phase of intense virus budding, i.e. host cells are in 
a late state of infection. However, the dynamics of the ceU 
proteome along the whole infection cycle may give important clues 
how virus components modulate the ceU proteome in the difierent 
phases of infection and, thereby, on early ceUular hubs as potential 
targets for inhibition of virus infection. 

Here, we present a quantitative study on the influenza A / 
Puerto Rico/8/1934 (HlNl) HlNl infection triggered time 
dependent dynamic changes of the proteome of infected Madin- 
Darby canine kidney (MDCK) cells monitoring for the first time 
viral and cellular proteins simultaneously. We assessed its impact 
on the host cell proteome by mass spectrometry based quantitative 
proteomics using the SILAC approach [21-23]. 

Results/Discussion 

To analyze dynamic processes along the whole viral infection 
cycle we generated quantitative data from the viral and the host 
cell proteome of influenza virus infected (MDCK) cells after 
various time points post infection (p.i.) by mass spectrometry 
analyses of SILAC samples. MDCK cells are of epithelial origin 
and serve as a model cell system for human epithelial cells since 
this cell line was used extensively to study influenza A virus [24] . 
The experimental procedure is illustrated in Figure la. In brief, 
MDCK cells were maintained in SILAC media containing 
variants of stable isotope labeled essential amino acids. The 
incorporation of these amino acids into newly synthesized proteins 
results in "light", "medium-heavy" and "heavy" cell populations 
[25]. Cells were infected with partially purified virus at an 
empirically determined amount of virus sufiicient to infect nearly 
all cells. HA titration of the stock suggested this was equivalent to 
an MOI of 100. Infection efficiency was ~98% as indicated by 
staining with antinucleoprotein antibodies (see SI and Figure SI). 
After infection the cells were harvested at various time points up to 
12 hrs p.i., digested by trypsin, and analyzed by mass spectrom- 
etry. This period covers also the onset of virus budding from 
infected MDCK cells, which occurs already at about 10 hrs p.i. 
[26-28]. The difierence in mass of the "hght", "medium-heavy" 



and "heavy" variants allowed the identification of the distinct 
fractions at once and thus, the direct comparison of the obtained 
data sets. The whole experiment was independendy replicated 
showing a good correlation and demonstrating the reproducibility 
of our assay (Figure S2). We identified 1311 MDCK cell proteins 
(Figure lb) and most of influenza A virus proteins (Table SI). In 
the following, we discuss the results for the virus proteins and for 
the MDCK cell proteins, as well as consequences for cellular 
regulatory and metabolic networks. 

Kinetics of virus protein expression 

The dynamics of the virus proteome were recorded up to 12 hrs 
p.i. and are presented in Figure 2. A few of proteins were not 
detected, for example M2, an alternati\ (; splicing variant of 
segment 7 [29]. As expected, the abundances of the viral proteins 
showed significant alterations throughout the time courses of the 
experiment. To provide independent evidence for the time course 
of virus expression measured by SILAC, we have recorded a 
immunofluorescence time course of NP expression (see Figure SI). 
NP expression becomes visible between 2 and 4 hrs p.i. with 
strong accumulation in the nucleus. Nuclear export was rarely 
observed at 6 hrs p.i., but became clearly visible at later time 
points. Expression reached a plateau around 10 hrs p.i. as also 
observed by the SILAC approach. 

Protein content fold changes are plotted for S4 hrs p.i. setting 
the amount at 4 hrs p.i. arbitrarily to 0. Note that most virus 
proteins exhibited an apparent decline in abundance between 0 
and 4 hrs p.i. which is due to our experimental set-up: The viral 
particles used for infection of the three differentially labeled cell 
populations were light. Combining samples is expected to cause a 
specific increase in light viral proteins with respect to the others 
and consequently an overestimation of viral protein abundance at 
the early phase of infection. Furthermore, this fraction of light 
virus proteins does not correspond to newly synthesized proteins. 
As a consequence, we have omitted plotting virus protein content 
for the early phase of infection. Theoretically, NSl should not 
decline since this protein is not part of incoming virus particles. 
Consistently, we observed that the decline of the protein is 
considerably lower than for the most other virus proteins (Table 
SI). The minor decline we detected might be due to measurement 
errors. A similar low decline was found for NEP which is only at a 
very low copy number present in virions. 

For >4 hrs p.i. a strong synthesis of all viral proteins was 
observed. Data for SlO hrs p.i. indicate that protein synthesis has 
reached a plateau. However, based on protein abundance two 
groups can be identified. Protein expression of HA, M 1 , NP, NSl, 
and - although somewhat lower - NA is much steeper (group 1) in 
comparison to NEP, PA, PBl, and PB2 (group 2). As outiined 
below the difierence between viral protein levels is in line with the 
function and the abundance of proteins in genesis and formation 
of new virions. For example, HA, Ml, NA belonging to group 1 
are the most abundant protein components of a virion. 

NSl shows a strong and continuous increase in abundance in 
the early phase of influenza A protein translation (Figure 2). This 
may ensure its efficiency in antagonizing the antiviral response of 
the host cell. As a non-structural protein it has been shown to 
influence mRNA splicing [30], mRNA nuclear export [3], mRNA 
adenylation [31] and translation [32,33]. In particular, suppres- 
sion of human mRNA nuclear export favors export of viral mRNA 
which is an important early step for synthesis of viral proteins. 
Furthermore, previous studies surmise that NSl impedes cellular 
antiviral defense mechanisms [6,34] by binding to a wide range of 
cellular molecules [3] . In fact, recent findings by Davey et al. 20 1 1 
[35] report about a mechanism that allows viruses to interact with 
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Figure 1. a. Outline of tKie experimental setup (For details see Material and Methods.], b, Proteomic phenotyping of the influenza 
A/PR/8 infected MDCK cell proteome using GO annotations. Quantiles of the quantification histogram are indicated at the top of the 

heatmap. Each quantile was separately analyzed for gene ontology pathways and clustered for the z transformed p values. The most prominent 
representatives of (Table S2) -represented biological processes of each quantile were selected and annotated in the right panel. 
doi:1 0.1 371/journal.pone.0094257.g001 



cellular proteins to hijack host resources. For NSl a short linear 
motif (SLiM) was identified suggesting this viral protein to mimic 
ceUular signal sequences and interact with cellular proteins. NSl is 
known to inhibit, to bind, and to inactivate the double stranded 
RNA dependent protein kinase (PKR) which is an important 
regulator of the anti-viral innate immune response in mammalian 
cells [36,37]. However, only recently it was found that NP of 
influenza A is also involved in inhibition of this enzyme but 
through a different mechanism [38]. 

Proteins of group 2 do not belong to components of the virus 
envelope. At this stage, we can only hypothesize why they are not 
present in high copy numbers in the initial phase of virus 
replication. For example, NEP together with Ml and cellular 
transport factors escorts newly formed vRNPs through the nuclear 
pore. It has been reported that NEP may shuttle back to the 
nucleus to catalyze export of further RNP complexes [39] which 
may provide an explanation for lower copy number of the protein. 
It is known that the polymerase complex factors PA, PBl, and PB2 
are expressed at rather low levels in comparison to other virus 
proteins [40]. Initial transcription cycles can be performed by the 
polymerase molecules coming with the infecting viral particles. 
After nuclear import, new polymerase subunits form trimers 
enabling further transcription and replication in particular at later 
stages of the infection [41]. Furthermore, at later time points, the 
newly synthesized polymerases will be complexed in vRNPs. 



The strong synthesis of NP and Ml could be related to 
resourcing of these proteins for formation of vRNPs, and, in case 
of Ml, also for formation of the budding site. NP is required in 
large amounts to form the vRNP complexes in concert with the 
viral antigenomic vRNA molecules [42,43]. Moreover, the level of 
NP is assumed to hold the key position for controlling the switch 
from transcription to rephcation mode [44]. This is supported by 
findings from a recent virus-host interaction screen [1 1] suggesting 
that NP might be involved in cellular signaling pathways and 
interfere with host cell responses, similar to the well described NS 1 
protein. Ml binding to those forming vRNPs is mandatory for 
subsequent nuclear export of vRNPs [45] . Only after formation of 
vRNP-Ml complexes, nuclear export becomes possible. This 
could rationalize the delayed steep increase of synthesis of NEP 
which is important for this export step [39,46] . 

Similar to M 1 , HA is also synthesized with fast progression after 
infection. As HA and Ml are transported to the plasma membrane 
of the host cell from which new virions bud, their steep increasing 
synthesis strongly supports the view that both are the key players of 
formation of the budding site. One reasonable scenario is that HA 
enriches locally in raft domains providing the docking site for other 
viral components. Several studies have suggested that the 
cytoplasmic tail of HA is recognized by Ml [47,48]. In 
comparison, NA synthesis is less pronounced. This could explain 
the lower number of copies in the envelope. It may also suggest 
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Figure 2. Time profile of viral proteins. Logarithmic presentation 
of the protein abundance fold change. For details see Material and 
Methods. 
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that a high abundance of this protein is not essential for the genesis 
of the budding site. Nevertheless, it seems to be of relevance in 
apical sorting of viral proteins [49] . 

In general, our data hint at a correlation between the 
abundance increase of viral proteins during the infection 
progression and their previously determined copy number in the 
viral particle. For example, the relative proteins levels of Ml, HA, 
and NA at 12 hrs p.i. (Figure 2) are in good agreement with the 
molecular composition of progeny viral particles [50] . 

Analysis of host cell protein dynamics upon viral 
infection 

In parallel, we measured and analyzed the host cell response up 
to 12 hrs p.i. to cover the dynamic changes of the host proteome 
composition along the time course from virus entry till release of 
progeny viruses. As a consequence of the integration of newly 
synthesized viral proteins into the data analysis which are not 
present at the starting point of the experiment, i.e. at very early 
stages of the infection, we quantified changes of the MDCK 
protein abundances relative to 8 hrs p.i., focusing on the 
identification of large-scale changes (Figure 3). To structure the 
huge data set (1311 identified proteins), we performed a fuzzy c- 
means clustering of the time profiles [51] and uncovered co- 
regulated sets of proteins. Gene ontology (GO) [52] and pathway 
(KEGG) [53] (Table S3) enrichment tests were used to relate 
protein sets to biological functions. Generally, proteins could be 
clustered according to the time dependent patterns of their 
abundances being either transiently or continuously up- or down- 
regulated. However, we did not detect an overall breakdown in 
protein expression as was suggested as consequence of the cap- 



snatching mechanism necessary for viral transcription initiation 
[54] . It is more likely that the virus is able to specifically control the 
cellular gene expression, since some proteins with antiviral 
function or those with functions less important for viral replication 
were suppressed and proteins that might support the production of 
progeny viral particles were increasing in abundance (see below). 

Although the expression profile is affected by virus infection, the 
overall protein expression is not changed much in the time course 
of our experiments. However, this does not preclude a redistri- 
bution of energy provided for different cell processes upon 
infection as indicated by the differential regulation of protein 
synthesis [28]. It must also be noted that our experimental set-up 
quantifies overall changes in protein levels that are affected by 
both synthesis and degradation. To directly quantify changes in 
protein production a pulsed SILAC (pSILAC) approach would be 
more appropriate [25]. 

Characterization of host cell protein dynamics upon viral 
infection 

With infection progression a downward trend for abundance of 
many host cell proteins was observed. Very likely, although 
ribosomal proteins are in general highly abundant and are 
expressed constitutively [55], accessibility of ribosomes for 
translation of cellular proteins becomes restricted due to enhanced 
translation of virus proteins. We surmise, eventually this wiU push 
the cell into an unsustainable state until the cellular biosynthesis 
collapses at the end of the replication cycle leading to the 
activation of apoptosis pathways. However, this warrants exper- 
imental validation. 

GO terms referring essentially to metabolic processes and 
protein synthesis in clusters 1 to 3 (Figure 3) are enriched in the 
early phase of infection. Anabolic and catabolic processes 
including glycolysis (Figure 3, clusters 1 and 2) show a strong 
increase up to 4 hrs p.i. followed by a temporary decrease 
approximately at 8 hrs p.i. in the SILAC profile (For additional 
information please refer to Figure S2.). A similar pattern was 
revealed for expression of proteins that relate to the "host anti- 
viral response" with GO terms such as immune response, 
histamine secretion and apoptosis (Figure 3, cluster 2). We 
observed a temporary down-regulation of this class of proteins 
approximately at 8 hrs p.i. This could reflect a virus induced 
antagonistic "host anti-viral response" activity as it is followed by 
the onset of production and release of progeny virions. In any case, 
we observed a strong reactivation of protein synthesis in clusters 1 
and 2 along with the onset and increased production of progeny 
virions (see also above). Remarkably, this increase is accompanied 
by reduction of ribosome biogenesis, rRNA processing, and 
translational elongation (cluster 6) indicating enhanced virus 
protein expression at the expense of expression of cellular proteins. 
Additionally, among these proteins characterized by a strong and 
constant increase in abundance were those connected to viral 
functions and replication cycle: GO terms enriched in cluster 2 
(Figure 3) include "localization and transport", "protein import 
into nucleus", "translational initiation", "actin filament organiza- 
tion" (actiii is involved in the formation of budding zones [56]) and 
"COPI vesicles" (required for post-translational modification of 
the viral surface proteins and their transport to the plasma 
membrane [57]). Another significant term was "response to heat" 
which has been found to be important in vRNA transcription [58]. 
Surprisingly, proteins that are constantly down regulated are to a 
large part related to DNA organization (Figure 3, cluster 4). 

We observed examples for differential regulation also within 
metabolism. Expression of proteins responsible for metabolic 
pathways with activation of glycolysis (Figure 3, cluster 2) behaves 



PLOS ONE I www.plosone.org 



4 



April 2014 I Volume 9 | Issue 4 | e94257 



Protein Profile of Influenza Virus Infected Cells 



Cluster dynamics 



Enriched GO terms 




oxoacidic metabolic process Q 
cellular catabolic process 
tRNA metabolic process 



glycolysis apoptosis Q 

localization vesicle transport 

ubiquitination histamine secretion 
proteolysis 



translational initiation © 
DNA dependent DNA replication 
reproduction 



TCA cycle 


o 


DNA packaging 




ion transport 





TCA cycle O 
oxoacidic metabolic process 



translational elongation Q 
rRNA processing 
ribosome biogenesis 
microtubule based process 



Oh 4h 8h lOh 12h 
Time (hrs) 



Figure 3. Functional clustering of all estimated cellular proteins. All normalized time profiles were clustered by a fuzzy clustering algorithm 
to find modules of co-regulated proteins. Enrichment tests for gene ontology terms on each cluster were performed for all proteins with a 
membership value >0.5 (n= number in brackets). The most significant terms are represented on the right panel (see also Table S3). 
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differently to that of proteins assigned to cellular respiration 
(Figure 3, clusters 4 and 5). Notably, the shift from respiration to 
glycolysis can also be found in proliferating cells or tumor cells 
[59-62] . Whether this dynamic energy compensation mechanism 
is actively controlled by the virus or the host cell itself remains to 
be elucidated. 

To get deeper insights into the diverse effects of the influenza A 
virus infection on the host cell proteome we employed proteomic 
phenotyping [63]. This method characterizes the differences 
between two proteome samples based on enrichment of classifi- 
cation terms. Profiling was performed using GO (Figure lb) and 
KEGG [64] (Figure 4) pathway annotations showing coherent 
results for the cellular proteome at 10 hrs p.i. where infection is 
already established but still below the maximum production of 
new virions. 

GO enrichment analysis revealed a strong increase of protein 
modules relevant (i) for viral replication itself, such as "proton 
transport", "RNA metabolism", "mRNA processing", "ATP 
synthesis", and (ii) for host cell response to infection related 
processes such as "cell death". The enhanced expression of 
proteins assigned to "protein localization" implicated an invasive 
response of the host cell towards viral infection 10 hrs p.i. This was 



also reflected by the inverse regulation of the ATP metabolism, 
which was significantly enriched in the up-regulated genes, but 
also represented in the down-regulated quantile. 

The proteomic phenotyping again revealed the switch from 
respiration to glycolysis. In the KEGG analysis (Figure 4) the 
"pentose phosphate shunt" and "ketone body formation" were 
enriched in the top 10% up-regulated genes, whereas for the 
"TCA cycle" higher z-scores were reached in the 10-50% range. 
On closer inspection, abundance of many of the host's glycolytic 
enzymes and those of the pentose phosphate shunt increased, 
while enzymes of the TCA cycle dropped slightly in concentration 
(Figure S3). This gradual shift from aerobic to anaerobic ATP 
production is in agreement with metabolic data gathered by Ritter 
et al. 2010 [28]. They reported elevated glycolysis rates and lower 
respiration after 8 hrs p.i. These changes are probably caused by 
the onset of apoptotic events in response to infection because the 
ceUular morphology is already changing at this phase of infection 
(data not shown). 

To include previous knowledge on host virus interaction in our 
analysis, we used gene interaction networks compiled by 
Watanabe et al. 2010 [65]. They combined data from various 
studies [7,8,1 1,66] to create networks of host virus interaction and 
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Figure 4. Proteomic phenotyping of the influenza A/PR/S infected MDCK cell proteome at 10 hrs p.i. using KEGG annotations. 
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to identify essential genes for virus replication. By combining these 
with our time-resolved data we were able to create a scheme on 
the viral influence on the host cell proteome. In our study we 
identified 26 out of 128 host proteins found in at least two previous 
screening studies and 80 proteins identified by interaction studies. 
Now, we can foUow the cellular changes taking place in each step 
of the replication cycle (Figure 5): (i) We detected an increase of 
abvmdance of the proton transporter ATP6V1E1 as well as of 
three other subunits of this transporter. ATP6 is essential for viral 
endocytosis [67,68] and acts in concert with HSPA8 (heat shock 



protein) mediating the uncoating of the clathrin coated vesicle 
which is important for endosome maturation and thus, the release 
of the viral genome from the endosome [69]. Notably, HSPA8 
showed an adequate qualitative behavior as ATP6. (ii) In the 
subsequent step the import of vRNPs into the nucleus is facilitated 
by the nuclear pore complex. Karyopherin (3 1 (importin) which 
was classified as essential for this active transport [70-72] and two 
additional subunits proved to interact with viral proteins were 
measured in our study and included in the combined network 
(Figure S4). They all showed similar activation dynamics during 
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For an efficient viral mRNA processing, a fine-tuned interaction 
between host and virus is crucial. The influenza virus pre-mRNA 
is processed using the cellular sphcing machinery and, in doing so, 
the detection by anti-viral mechanisms of the host is circumvented. 
We detected 14 proteins, among them RBMX, RBM39, and 
RBM25, that are related to mRNA processing and interact with 
viral proteins [65] . (iv) The translation of viral proteins is strongly 
dependent on the cellular biosynthesis machinery. We identified 
26 of the 42 ribosomal proteins that were reported to interact with 
vRNP [73]. Contrary to the general trend of reduced protein 
abundance, seven out of these proteins increased with progression 
of infection. Among all activated proteins which we have 
identified, three ribosomal proteins have been reported to exhibit 
extra-ribosomal functions in p53 activation and inhibition of 
proliferation [74]. This is in accordance with the observed shift 
from respiration to glycolysis and indicates an overall change in 
the cellular status. Further, an activation of translation initiation 
factors was detected, such as EIF2, EIF3 and EIF4, which are 
essential for efficient viral replication [75]. The positive influence 



of virus infection on the expression of EIF4A2 known to enhance 
the cap recognition at ribosomes supports the previous observation 
on the enormous efforts towards viral translation [76]. 

In conclusion, considering the impact of the infection on the 
cellular protein expression profile and the consequences for 
various cellular networks one can hypothesize that the virus not 
simply hijacks existing resources - it modifies the whole cell status 
for its own purpose. It becomes evident that both, the response of 
the host cell and the emergence of newly synthesized virus 
components, resemble a dynamic system which would be difficult 
to overview and to understand by investigating its numerous parts 
separately. As we have shown here, our and complementary 
approaches enable to follow the dynamics of the viral and cellular 
proteomes quantitatively and the consequences on cellular 
networks along the entire virus infection cycle. This proteomic 
approach on a systems level paves the way for new strategies to 
combat viral infection. For example, extending this holistic 
approach by comparison of the proteome pattern of infected host 
cells between different influenza virus strains and correlation with 
the respective infectivity should highlight main checkpoints of the 
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Figure 5. Schematic depiction of the steps of the viral life cycle and the associated host, proteins. Proteins were selected from the 
overlap between proteins identified by our approach and genes classified as essential for viral reproduction or interacting with viral proteins by 
Watanabe et al. 2010 [68]. These proteins were grouped according 453 to their function and put into context by the schematic graphics of the 
infection cycle. 

doi:1 0.1 371/journal.pone.0094257.g005 



PLOS ONE I www.plosone.org 



7 



April 2014 I Volume 9 | Issue 4 | e94257 



Protein Profile of Influenza Virus Infected Cells 



Table 1. Scheme of the influenza A/PR/8 infected MDCK cell 
sample preparation. 





Variant 


Harvesting time 
point (hrs) 


MDCK cells in light SILAC medium 


0, 10 post infection 


MDCK cells in medium heavy SILAC medium 


4, 12 post infection 


MDCK cells in heavy SILAC medium 


8 post infection 


doi:10.1371/journal.pone.0094257.t001 



host cell essential for a permissive infection. Subsequent studies 
will unravel whether those checkpoints may serve as targets for 
antiviral drugs. 

Material and Methods 

Cell Culture, virus, and sample preparation 

Mammalian MDCK cells were cultivated in DMEM supple- 
mented with 10% FCS in a humid incubator at 37°C and 5% 
COj. Due to the metabolic incorporation of stable isotopic 
variants of amino acids, the applied FCS had to be dialysed to 
eliminate the natural amino acids. SILAC media were prepared as 
described previously [77] "Heavy" and "medium-heavy" SILAC 
media were prepared by adding 84 mg/1 '^C6''^N4 L-arginine plus 
146 mg/1 '''C6''^N2 L-lysine and by adding 84 mg/1 '"^Ce-L- 
arginine plus 146 mg/1 D4- L-lysine, respectively. Labeled amino 
acids were obtained from Sigma Isotec ( Cg-L-arginine, C^ ' N4 
L-arginine and ''^C6''^N2 L-lysine) and Cambridge Isotope 
Laboratories (D4-L-lysine). To prepare "light" SILAC medium, 
the corresponding non-labeled amino acids (Sigma) were added. 
MDCK cells were split every 3rd to 4th day in a 1:10 ratio. 
Resulting from cultivation in the respective SILAC medium over a 
period of at least eight passages MDCK cells incorporated 
uniformly the isotopic amino acids in their proteome. 

Cells subjected to virus infection were seeded in 6 well plates at 
80% confluency in DMEM overnight. Influenza A virus /Puerto 
Rico/ 8/1934 (HlNl) was produced in 10-day-old chicken 
embryos and purified as described previously [78]. The allantoic 
fluid was collected, and cell debris was removed by a low speed 
spin (880 x^, 30 min). The virus was pelleted by spinning the 
allantoic fluid at 95,000 for 90 min. The pellet was resuspended 
in PBS and homogenized with a Teflon-coated homogenizer. The 
HA titer was determined by hemagglutination of red blood cells 
(RBCs). One hemagglutinating vmit corresponds to a titer of 
1,6x10 6 plaque forming units. 

In all indicated experiments infection with influenza A virus/ 
Puerto Rico/ 8/1934 (HlNl) was performed at a multiplicity of 
infection (M.O.I.) of 100. To this end, the virus solution was added 
directly into 1 ml DMEM and added over the cells. After one hour 
incubation at 37°C to allow the virus to attach to the cell surface, 
the supernatant was replaced by 1 ml fresh DMEM and 
incubation was continued. At various time points post infection 
(indicated in table 1) cells were harvested by trypsination and 
pooled according to the scheme in table 2. After cell pellets were 
stored at — 20°C for further procedure. 

Cell lysis and preparation for in-solution digestion 

Cells were lysed in appropriate amovmts of denaturation buffer 
(6 M urea/2 M thiourea in 10 mM HEPES, pH 8.0) for 20 min 
on ice. The lysates were cleared by centrifugation for 10 min 
(14,000 rpm at 4°C) and transferred to fresh tubes. Protein 



Table 2. Scheme of samples subjected to mass spectrometry. 





Sample ID 


Composition 


Sample 1: 


0 p.i. (light) +4 hrs p.i. (medium heavy) +8 hrs p.i. (heavy) 


Sample II: 


10 hrs p.i. (light) +12 hrs p.i. (medium heavy) +8 hrs p.i. (heavy) 



doi:l 0.1 371 /journal.pone.0094257.t002 



samples were reduced for 30 min at RT in 10 mM dithiothreitol 
solution followed by alkylation for 20 min by 55 mM iodaceta- 
mide in the dark at RT. The endoproteinase LysC (Wako) was 
added (proteinienzyme ratio of 50:1) and incubated for 4 hrs at 
room temperature. After dilution of the sample with four times 
digestion buffer (50 mM ammonium bi-carbonate (NH4HCO3) in 
water, pH 8.0), sequence grade modified trypsin (Promega; 
protein: enzyme ratio of 50:1) was used for overnight digestion. 
Trypsin and Lys-C activity was quenched by adding trifluoroacetic 
acid to adjust the pH to <2, and peptides were extracted and 
desalted using StageTips [79] . 

HPLC and mass spectrometry 

Reversed-phase liquid chromatography (rpHPLC) was done on 
an Eksigent NanoLC - ID Plus system using self-made fritless C18 
microcolumns [80] (75 |im ID packed with ReproSU-Pur C18-AQ^ 
3-|im resin. Dr. Maisch GmbH) connected on-line to the 
electrospray ion source (Proxeon) of a LTQ;Orbitrap Velos mass 
spectrometer (Thermo Fisher). Peptide samples were picked up by 
the autosampler and loaded onto the column with a flow rate of 
250 nl/min. Subsequent sample elution was performed at a flow 
rate of 200 nl/min with a 10 to 60% acetonitrile gradient over 
6 hrs in 0.5% acetic acid for online MS analysis. The LTQ_- 
Orbitrap Velos instrument was operated in the data dependent 
mode (DDA) with a full scan in the Orbitrap followed by up to 20 
consecutive MS/MS scans in the LTQ. Precursor ion scans (m/ z 
300-1700) were acquired in the Orbitrap part of the instrument 
(resolution R = 60,000; target value of 1 x lO''), while in parallel the 
20 most intense ions were isolated (target value of 3,000; 
monoisotopic precursor selection enabled) and fragmented in the 
LTQ, part of the instrument by collision induced dissociation 
(CID; normalized collision energy 35%; wideband activation 
enabled). Ions with an unassigned charge state and singly charged 
ions were rejected. Former target ions selected for MS/MS were 
dynamically excluded for 60 s. Total cycle time for one fuU scan 
plus up to 20 MS/MS scans was approximately 2 s. 

Processing of mass spectrometry data 

Identification and quantification of proteins were performed 
with the MaxQuant software package [81]. The software 
automatically extracts, re-calibrates and quantifies isotope clusters 
and SILAC triplets in the raw data files (heavy labels: ArglO and 
Lys8; medium-heavy labels: Arg6 and Lys4; maximum of three 
labeled amino acids per peptide; polymer detection enabled; top 

6 MS/MS peaks per 100 Da). The generated peak lists were 
searched against a target-decoy database of forward and reversed 
proteins [82] by the built-in search engine Andromeda [83]. The 
search was done using an Influenza A virus /Puerto Rico/8/1934 
(HlNl) database (obtained from www.biomart.org in July, 2010) 
and a domestic doc {Canis familaris] [84] database. FuU tryptic 
specificity was required and a maximum of two missed cleavages 
as well as a mass tolerance of 0.5 Da for fragment ions was 
allowed. The initial mass accuracy cut-off on the precursor ion was 

7 ppm but subsequently narrowed down by filtering based on hits 
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to reversed peptides in tlie target-decoy database. Oxidation of 
methionine and acetylation of the protein N-terminus were set as 
variable modifications, carbamidomethylation of cysteines used as 
a fixed modification. For filtering of peptide identifications, 
assembly of proteins and re-quantification the following param- 
eters were used in MaxQuant: A minimum peptide length of 6 
amino acids was required and the false discovery rates (FDR) were 
estimated based on matches to reversed sequences in the 
concatenated target-decoy database. A maximum false discovery 
rate of 1 % at both the peptide and the protein level was permitted. 
Peptides were assigned to protein groups (that is, a cluster of a base 
protein plus additional proteins matching to a subset of the same 
peptides). Protein groups containing matches to proteins from the 
reversed database or contaminants were discarded. Protein ratios 
(M/L, H/L, H/M) were calculated as the median of all peptide 
ratios for a specific protein. Only proteins with at least three 
peptide ratios at each measuring time point were considered. 

Data analysis 

For the clustering we standardized the fold-changes of the time 
course to mean zero and standard deviation of one. A soft 
clustering algorithm calk'd fuzzy-c-means from the Mfuzz [51] 
package in Bioconductor [8,5] was used. We chose this method, 
because it is very robust against noise and provides an easy way for 
posterior filtering of important genes in the different clusters. The 
method appoints a membership value between zero and one for 
each protein to each cluster. The parameter m determines the 
influence of noise on the establishment of cluster centers. Higher 
values of m decrease the influence of poorly classified data points 
on the cluster centers. Applying an iterative approach the number 
of cluster centers c = 6 and the parameter m = 1.8 were defined. 

Identified proteins were assigned to their biological process 
using gene ontology (GO) [52] enrichment analysis. Enrichment 
for biological process (BP) was tested using a hypergeometric test 
from the GoStats [86] package in Bioconductor for the single 
clusters with the complete set of measured genes as a background 
distribution. This test automatically corrects for the bias resulting 
from the tree structure of the ontology. The analysis was as well 
performed for the KEGG [53] pathway annotations in a similar 
fashion. To characterize the proteomic changes provoked by the 
influenza A virus infection a proteomic phenot)'ping for GO terms 
at the time point 10 hrs p.i. was performed as previously described 
[87]. We divided the skewed distribution of measured log2 fold 
changes into four quantiles (10%, 50%, 90%, and 100%) and did a 
GO enrichment test for each of the quantiles with all measured 
genes as background. For all GO terms that appeared in one of the 
quantiles with a p-value <0.05 we did a transformation by - 
loglO(p-value) and computed a z-score by (x-mean(x))/sd(x). 
Colors in the heatmaps encode the value of these z-values. 
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